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ABSTRACT 



The effort here was concerned with initiating a long 
term experimental study to develop the method of laser 
light scattering as a plasma diagnostic technique. This 
was done with a ruby laser and a grating spectometer cou- 
pled with sensitive detectors, at the steady state plasma 
facility at the Naval Postgraduate School. Although at- 
tempts to detect the extremely low levels of scattered 
light failed, excellent progress was made in determining 
existing apparatus limitations and in recognizing which 
system parameter may be adjusted in order to detect the 
scattered light; this will enable the future experimenter 
to progress with a minimum of overlap. 
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I. 



INTRODUCTION 



This thesis is intended as a detailed report of the 
beginning of a long-term experimental effort to develop 
Thomson scattering of laser beams as a plasma diagnostic 
technique. The technique employs the particle-like behav- 
ior of photons in scattering off plasma electrons. In 
general the background radiation from an energetic plasma 
dominates making it impossible to distinguish the scatter- 
ed light. This problem may be somewhat alleviated by the 
use of an extremely intense, sharply tuned source in con- 
junction with narrow band detectors. The high power ruby 
laser is such a source and a grating spectrograph-photo 
multiplier tube arrangement is such a detector. Even with 
these improvements, a manifold of problems and complica- 
tions yet exist, not the least of which is careful planning 
of the geometry of the experiments and associated alignment 
of apparatus. Previous experiments have shown that laser 

scattering from electron beams of energy approximately 2 

9 -3 

kev and density approximately 5 x 10 cm yielded scatter- 

— 18 

ing intensities of only about 10 of the laser output. 
Though electron densities in the plasma under study are 
significantly higher, this gives some idea as to the prob- 
lems presented in discerning the scattered light. 
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It will be shown that the spectral distribution of the 
scattered light is in part determined by a dimensionless 
coherence parameter proportional to the ratio of incident 
light wavelength to the Debye length; 

ff oC i- 
o 

when ol is small, that is, when the wavelength of incident 
light is small compared with the Debye length the cross 
section is the appropriate Thomson cross section and the 
shape of the scattered spectrum is a Gaussian; the width 
is the Doppler width associated with electron motion and 
as such is a measure of electron temperature. 

When a is large (ot » l) the spectrum has a central 
line of the incident light frequency, , between two sat- 
ellite peaks at frequencies 03 + U) where o,' is the plasma 

op p 

frequency. In this case the Doppler broadening of the cen- 
tral peak is approximately proportional to the ion motion 
and gives information concerning ion temperatures. 

The research begun in the plasma laboratory of the 
Naval Postgraduate School has been under the cognizance of 
Professor F. Schwirzke and is financed under Naval Air Sys- 
tem Command Research Grant. The accomplishments made thus 
far and described herein will enable to the future experi- 
menter to progress with a minimum of overlap. The major 
obstacle yet to be overcome, and indeed the main problem 
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in the development of this diagnostic method as a whole, is 
detecting the scattered light and distinguishing it from 
natural background radiation. Once this is done, the tech- 
nique may be rapidly and effectively developed. 

The diagnostic technique once developed, it may be 
employed to investigate turbulence in a laser produced 
plasma. 
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II. METHOD AND THEORY 



A. ELECTROMAGNETIC THEORY OF LIGHT SCATTERING 

Following directly from the equation of motion the am- 
plitude of an electron oscillation caused by the electric 
field E 1 of the incident electromagnetic wave of frequency 

W o is * 



S(t) 



e E* ( t ) eX 2 E 1 (t) 

2 2 2 
m 60 me 

o 



(1) 



Multiplying the displacement by the electron charge e 
the induced dipole moment is obtained; 



p = - es = 



-e 2 X 2 E ! 

4ff 2 mc 2 



€E' 



( 2 ) 



For the laser light scattering to be dealt with here, 
is a plane monochromatic wave cut off exponentially at 



time T = 1 / y 


, the pulsetime constant. 




That is, 






E' (r, t) 


= E' exp[i[k'»_r + (60 c 4 iy) t] } 


(3) 


P(r, t) 


= eE'exp[i[k’ «r 4 ( 60 o 4 iy) t] } 


(4) 



It is this time dependent Dipole moment (4) which produces 
the scattered radiation. 
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electron 



The electric field Ej"(r,t) scattered by the j ^ 
at position r^ detected at ft, from the Hertz formula for 
dipole radiation, is 



E_. " (ft, t) 



exp[i(k' -j:_. + k'R^) 



(5) 




( 6 ) 



where is the distance between scattering center and de- 



tec tor 



f^j(Rj,t) is the amplitude of the scattered field and is 



f ■ (R.,t)= 4 ” 9 [R. x (R. x €E’ ) ]exp[io> -y) t] 

J 3 R .V ^ 3 o 

3 



(7) 



that is, € = 



2.2 

e A 

2 2 
477 me 



) 



fj(Rj/t) = - (-^ 2 ') x(a.x E')]exp(iw t-yt) (8) 



R . — 1 —j 

me “j 



where a . = j — is the unit vector in the R. direction 

1 I Rj | 3 

r^ is measured from a reference point, 0, in the scattering 
volume. If the dimensions of the scattering volume are 
small compared with ft-, then R^ may be replaced by ft ^ , and 
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f_j = f for all electrons in the plasma volume. When deal- 
ing with the phase, however, more care must be taken in 
approximations. The following substitution may be made for 
Lj K< • 



R. = </(R-r ,) 2 = y ( R 2 -2R* r . +r . 2> « + ... ( 9 ) 

— j j — —j j st v ' 

Writing the scattered wave vector k" as having approximately 
the magnitude of the incident wave vector k' and direction 
that is. 




(5) may be rewritten as 

Ej (&, t) = j£(_0, t, X0) exp (ik ' fi. + ik*iTj) 
where k is the scattering vector. That is 

k = k"-k ' , k = (^ y^- sin y , k ' sa k" = y^ 



( 10 ) 



(ID 



( 12 ) 



The scattering angle 0 is between k' and k " . The direction 
is of course determined by the incident laser beam and by 
the position of the detector. 

In 90° scattering, 

T 

I 

I 




da 
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Furthermore X'« X" = X. The wavelength difference between 
X' and X"’is usually small and may be neglected, such that 
the symmetric of the scattered spectrum is undisturbed; if 
is the frequency of the incident light and is the fre- 
quency shift of the scattered light, I ( u^+CO) « I ( CO^-U!) if no 
drift velocities exist. 

The phase of the scattered light depends on the posi- 

ik- r . 

tion of the scattering center by the factor e . This 

becomes important when the fields scattered by a number of 
centers are superimposed at Slamming, over all electrons 

in the scattering center, one obtains 



E" (R, t) 



^Ej(R, t) = f(R, t)e lk ’ R 'le j 

j=l j=l 



( 13 ) 



The scattered intensity is equal to the time averaged 

I II I 3 Q 

Poynting vector | IC | ; Introducing the incident light 

intensity as 




one obtains 



( 14 ) 



1(0, X) 



IOO 



N 

CT ( 0 ) < | ^ e 

j=l 



i (k* r . ) 

— r 




( 15 ) 
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Where a (8 ) , Thomson's scattering factor, is equal in essence 
to | f 2 |; Introducing e‘ and e", the polarization vectors 
of the incident and scattered light, 

E ' = e'|E' | , E" = e" | E" | , ft-e" = 0 (16) 



o(B) 




(ee 1 ) 



4ff 

X 2 Si 



(17) 



for linearly polarized light; and for natural incidents 
light. 



' (T(0) 




(1 + cos 9) 
2ft 2 



(18) 



where the 1 in parentheses represents the contribution 
from natural incident light and cos 0 represents the con- 
tributed polarized parallel to the plane of scattering. It 
is easily seen how (17) complies with this. 

B. SPECTRAL INTENSITY DISTRIBUTION OF SCATTERED LIGHT 

For the most part, it is the electron density fluctu- 
ations in the plasma which are responsible for the scatter- 
ing of light; that is, according to Huygen's principle no 
scattering can occur in a strictly homogeneous medium. The 
electron density n (j:, t) is a function of position and time, 
and may be considered as a continuously varying quantity or 
as a superposition of delta functions as circumstances 
warrant . 
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Imagine the plasma confined to a large cubical volume 
3 

V = L and write the electron density as a spacial fourier 
1 

series . 



n (£'t) = ^ £ n^tje 1 ^ 

q - 



( 19 ) 



Where q is the wave vector satisfying the cyclic boundary 
conditions 



_ 2nhi , 

q -i “ h 



= 0, ±1, ±2 , h i = x, y, z 



( 20 ) 



The fourier coefficient n^ may be considered as the ampli- 
tude of the longitudinal plasma wave with wave vectors q 
and may be written 




n(ry t) exp (-iq.rjdr 
V 



for the microscopic representation and 



( 21 ) 



N 

n q = I exp(-iq-rj) 

j=l 



( 22 ) 



for the microscopic representation. 

Usually, longitudinal oscillations in plasmas have com- 
plicated time dependencies and dispersion relations; further- 
more a number of frequencies may be present. The probability 



1 

As done in Light Scattering and Drift Instability in 
a Plasma by Otto Theimer. ltc. 
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of a certain frequency occurring is related to the mean 

square of the amplitude n (to) of a wave with wave vector 

Si 

q and frequency 00 . A long lived, relatively undamped 
oscillation will have a sharp maximum in n c (oo) . A con- 
tinuous function of the frequency n (0>) is the Fourier 

Si 

transform of nq(t) with respect to time and may be written 
in integral form; 



n q (W) = J n q (t) exp(-iu>t)dt (23) 

But the actual density fluctuation is not by itself of pri- 
mary concern in scattering theory; rather, the amplitude 
of the polarization waves Pg(t) induced by the incident ra- 
diation is of interest. [Not to be confused with the dipole 
moments of equation (2)3 

P (t) = € n (t)E’ (t) (24) 

Si SL 

Where e is given by (7a) and is the polarization of the 
free electrons. Again, 12 (t) is the electric field of the 
incident laser pulse and has the time dependence 



E/ ( t) = E ' exp[ (iW Q -y) t] , t > 0 

e' (t) = 0 , t < 0 



(25) 



That is, the laser is fired at t = 0 and is essentially 
overafter T = ~ . Thus it is seen that the polarization 
waves have a frequency dependence (ci) c +a!) the amplitude of 



16 



which decreases exponentially with the laser pulse, as 
exp(-yt). That scattered light will of course have the 
same time dependence, and the scattering waves may be iden- 
tified as a frequency shift in the scattered spectrum. 

Thus the finite laser pulse may be attributed (at 

least mathematically) to the plasma waves and n (CO) may be 

Si 

written as 



n (co) = n (t) ex P C-(iw+ y) t]dt 
S J o - 



(26) 



Considering the total scattered intensity, the time aver- 
age value of the function n^ft) is taken 

n k ( t) I 2 exp(-2yt)dt 



< CIV t) l > ' ~ 2 



oo- 



exp (-2yt) dt 



00 

= 2y I k (t) \ exp (-2yt) dt 



(27) 



By Parseval 1 s Theorem from (26) and (27) it follows 

that 



<K(t)| 2 > = ? ,[_ ro K<«) I 2 dCO . (28) 



Clearly, additively combining the square of the amplitudes 
(ff) I i\(co) I gives the square of the amplitude n^(t)| , 



In Bronstein and La. Semgndjajew, Taschenbuch d. Mate 
matik,Verlag Harri Deutsch, Frankfurt/M. , 1965, va5 pp. 404-474 
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Thus each wave of wave vector k and frequency 0) contributes 
to the total scattered intensity by an amount given by 



I (co ,k) dco = 



CT ( 0 ) Cl ) I n k^ 



watts 1 
cm 2 -* 



( 29 ) 



And as already stated, each frequency shows up as a fre- 
quency shift in the scattered spectrum. 

So it has been established that each electron becomes 
Hertz oscillator, when induced to oscillate by the inci- 
dent light, and itself emits the light which is the scatter- 
ed radiation. If all electrons oscillate in phase with the 
incident light, the emitted (scattered) radiation will be 
coherent in time. If, however, the electrons have complete- 
ly arbitrary spacial distribution, coherence is destroyed. 
But electrostatic interaction prevents the electron distri- 
bution from ever being completely arbitrary, and it is from 
partial spacial coherence that so-called collective scat- 
tered radiation is observed. The degree or intensity of 
coherence is dependent on the ratio Ot between the wave- 
length and the characteristic Debye shielding length for 
the plasma. That is. 



a 

s 



X 

— x — : — & 

4n X sin 2 



( 



KTs 
47 t ne^ 




s 



i, e 



( 30 ) 



Where X is the wavelength of the incident light, X^ is the 
Debye length and 0 the scattering angle. For a « 1 scat- 
tering is incoherent; for a » 1 spacial coherence 
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determines the shape of the scattered spectrum. Perhaps a 
more physical and qualitative (if not mathematically rigor- 
ous) understanding may be had of this concept if one imag- 
ines the Debye spheres as integral units. Imagine the elec- 
trons as tied to the ions they shield. Within the Debye 
sphere, the electron distribution is arbitrary; that is, 
incident of wavelength less than the radius of the Debye 
sphere will oscillate electrons within that sphere in dif- 
ferent phases, whereas incident light of wavelength greater 
than the Debye sphere will oscillate the Debye spheres, and 
thus all electrons within them in unison, in coherence. 

The incoherent scattering spectrum is a direct illus- 
tration of the electron velocity distribution, as will now 
be shown. As the scattering under consideration is inco- 
herent, the effect of a single electron may be looked at. 

It has been seen that the electric field of light scattered 
fcll 

by the j electron, by equations (7) and (11) is propor- 
tional to the phase factor expCia^t + ik-r^]. Approximat- 
ing jTj by an expression valid for straight line electron 
paths , 

Lj ( fc ) = Lj (°) + Vj t , (31) 

the scattered field has the frequency 

0>" = CO + k»v. (32) 

o j 

and the scattered light is Doppler shifted by k*Vj = kv^ 

(k = k"-k'). The shift is two fold; that is the electron 
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is a moving source to the scattered light. Thus the total 
shift is k*v = k" • v - k'*v. F(v) now, is the normalized 
velocity distribution function with velocity v parallel to 
k. Corresponding to velocity v, the frequency is 40 = kv; 
the velocity interval dv correspond to frequency interval 
d40 = kdv. Thus, the number of electrons with velocity v 
is proportional to the scattered intensity with shift kv. 
That is. 



F e ( v) 



kl (kv,jj) 
I(k) 



(33) 



where v is the velocity component parallel to k, I(kv,k) is 
the scattered intensity integrated over all frequencies. 
Again, this is only valid for a « 1. 

The topic now turns to the more complicated theory of 
longitudinal plasma waves. The attempt here is to calculate 
the mean thermal value of the square of the space time de- 
pendent Fourier coefficients | n^(to)|^. The basic reference 

— 3 

for this calculation is the paper by Salpeter. The ampli- 
tude squares | n^ (40) | ^ will be solutions of the collisionless 
Boltzmam equation, i.e. the Vlassov equation (See Appendix A 
for derivation of the Boltzmann equation) . The collision- 
less Boltzmann equation. 



5f 

at 



+ v * V f 
r 




V f 

— v 



0 



(34a) 



3 

E. E. Salpeter, Phys. Rev 12, 1528 (1960). 
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is essentially a continuity equation for representative 
points in electron phase space. The increase in time in 
the number of electron points in phase space is caused by 
three factors: (1) The drift in points in position coordi- 
nates due to inherent velocity: (2) the drift in points 

in velocity coordinates due to electric and magnetic forces, 
and; (3) the spontaneous velocity change without position 
coordinate due to collisions. 

In the purpose of calculating | | the spacial 



Fourier coefficient f (t) of the distribution function 



f(r^,v, t) is introduced. (q is the general wave vector). 
Corresponding to equation (22) , 



the sum being taken over those particles with velocities 
between v+dv. As dv “* 0, the right hand side has one term 
if a particle of velocity v exists and no term if such a 
particle does not exist. Thus (34) represents a fine grain - 
ed distribution function so described because the density 
has but two values: zero in a vacuum and infinite at the 
position where a particle is represented by a mathematical 
point. Clearly, then, an integral over the velocity dis- 
tribution becomes identical with a sum over all the parti- 
cles : 




qz 



VUV ) 

expC-iq* r j ( t) ] , 

j 



(vdv) 



(34) 




(35) 



j=l 



21 



The number of terms in (34) is equal to the number of par- 
ticles with velocities in _v and _v + d_v. If the physical 
properties that are influenced only slightly by the micro- 
scopic details of the distribution function are calculated 
and if the course grained (continuous) distribution is in- 
dependent of spacial position, the number of terms in the 
sum is the most important aspect of (34) , and may be approx- 
imated by 

f qv (t) = n(r, t)F(v) = V -~ ~ - 2, n q ( t) exp (iq* r.) (36) 

q ” 

where F(v) is the normalized coarse grained (continuous) 
velocity distribution and n q (t) is the Fourier coefficient 
of the particle density given by (21) and (22) . As the no- 
tation in (36) implies, F ( v) is independent of time. This 
of course, is valid only if Af/F is small during the time 
l/Ao> where Aco is the practical spectral resolution. 

Without regard to the coarse or fine-grained structure 
of f (t), the distribution function may be represented for 
ions and electrons as a Fourier series; 



P g (r#v, t) = ^ Y f s qv (t) exp(iq.r) , s = i,e 



(37) 



Similarly, for the three terms in the Vlassov equation, 



9f s 1 V 9f sc<v( t ) 

1F = V A exp(iq-r) 



( 38 ) 
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and for the second term 





exp(iq*j:) 



(39) 



In the third term, only the electrostatic force e E (r_, t) is 

s 

considered where e g is the charge of the particle at r-E(r,t) 
is written as a Fourier series. 



E(r,t) = ^ (t) exp(iq.r) 



V L> =q' 

q - 



(40) 



For equation (36) is used, where F(v) is the only ve- 

locity dependent quantity. The third term then becomes. 



$ Sf e 
s s s 



m ’ 5v „2 ^ l ‘ ii q 

s - m s V q 2 



V r Sf (v) 

2,[E . exp(iq' .r>2n „exp(iq".r) ] 1- — (41) 



q S 2 



Combining terms in the double sum for which 



q + q" = q » q" = q - q' 



(42) 



9 df e 
s s s 



SF s ( v) 



m civ ..2 

s - ra s V q q 



nX' n s<q-q') eX P (i a-£>- — 57 

rr Cf — — — — 



(43) 



1 V 3F s 

= —2 l Z ex P (i S'^) * 3v 



m V 



a. - 



where 



y E , n 

i Is —a s 



^ — e / s . ii , »\ 

sq S q«~S s(q-q') 



( 44 ) 
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The "Random Phase" approximation is now introduced. The 
order of magnitude of n^ follows from (22) . Remember that 



Lim n = N » n - JW 
N-“ 



(45) 



n^ is the sum of N positive and negative term with values 

arbitrarily distributed from zero to one. The probability 

that the real and imaginary parts of n have the values 

£ and V respectively is 

£L SL 



U CT + ) 



TT/t » v _±_ r - 'an 

W ^q ,T7 q' 21TN SXp L 2N -1 



(46) 



and the mean amplitude square value is 



CO CO 



<1 n a | 2 > = i _ + ’ , a 2)WU a'’ , q ) = 2N (47) 



q q' 



(46| and (47) explain the importance of (45). E q that is, 
q = 0 has no dominating position for the electric field 
because it has a mean value in a neutral plasma of zero. 
Thus (44) has one dominating term only. 

Writing explicitly. 




e 

s 



|E K + ) E ,n , ,,! 

L— q L — q' s(q-q')J 

q’*q 



(48) 
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The "random phase" approximation assumes that the sum can 

4 

be neglected in comparison to the term E^ N. Adopting this 
approximation, equation (33) finally becomes 



^ Jexp(iq-r) [- 



c>f 



~sqv 

St 



+ 



ivqf + 
— _ sqv 



e n 
s s 



m 



E • 

— q 
s — 



SF (v) 
s — 

Sv 



]■ 



o 



(49) 



The Fourier representation is unambiguous, and thus (49) is 
valid not only for the sum of all general wave numbers q 
but also for every single q . Writing (49) explicitly for 



electrons and multiplying by the 'electron charge e, the 



equation for the Fourier coefficient cr corresponding 

e • £> v 

to charge density cr : 

J J ev 



S a 

eqv 



at 



+ 



iv* q 



a 

eqv_ 



2 

e n 

+ E . 

m — q 
e _ 



SF (v) 
9v 



0 



(50) 



An analytical expression for E may be obtained from 
Poisson's equation 

V • E = 4frcr = 4 ir(a +cr.) (51) 

— — e i' ' 

Where a is the charge density created by electrons and ions. 
As previously done, let 

CT(r,t) = ^ X a q (t)exp(iq-r) (52) 



4 

Though the infinite number of terms in the sum makes 
it by no means really clear that the approximation is al- 
lowable . 



25 



From (52) , (51) » and (40) 



I 



iq*E 
“ ~£ 




-4 fficr 



q 



q 



and q = k (as in equation (12) ) 



(53) 



Ber 



ekv 

Bt 



+ ikv c r 



ekv 



z 4rrn 

( e 



m 



ek 



Os 



BF e ( v) 
3v 



(54) 



Where v and E are the components of v and ]3 parallel to k 
(k is along the z axis) . The Fourier time transform of 
°ekv ; '' S neec ^ ec ^ i- n t ^ ie calculation of the spectral distribu- 
tion; this is written as in equation (20) 






CT ekv^ ^ ex P^“( ia)+ ^ ) t ^ dt 



(55) 



Integrating over all velocities. 



CT ekv (a,)dv S Q e k^ = J <^ ek ( t) exp[- (ico+y ) t]dt (56) 

— 00 o 



q e k v ( w ) satisfies partial differential equation obtained by 
equating the transform on the right and left of (54) . The 
following integral is needed 



09 Bcr 



exp[- (ico+y) t]dt = -c ekv (0) +(ico+y) q ekv (co) (57) 
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This is obtained by integrating by parts and using (55) . 



Notice that the contribution from the upper limit cr e ^ v ( co ) 
disappears due to exp(-yt). Using (57) in (54) gives 



, 4ffn dF (v) 

< W a ’ )(i “ +y+ikv)= C 'e*v (0)+ 1 



ek 



dv 



(58) 



or multiplying by i ~ 

^ 4ffn N dF (v) 



ekv 



ek 



dv 



q ekv ( “> = 



-co - kv + iy 



(59) 



where 



= Qek ( “> + °ik<"> 



(60) 



The charge density for the ions is obtained similarly, that 
is. 



lOr .. 
lkv 



^4lTn >. dF. (v) 

(0) -wf“) Q k ( ^ -kr 



‘J i kv ( “ ) = 



ek 



-co - kv +iy 



and the total charge density 



(61) 



q kv ( “ )=q ekv (t0)+q ikv ( “ ) = 



iCT kv<°>- 




-co - kv +ico 



s F ± (y)-| 

dv J 



(62) 

Integrating (62) over all velocities gives an expression for 
Q^(CO), which may then be introduced into (59) and that 
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result integrated over all v. The final result would be 
the amplitude necessary for scattering = - en^(W) . 

Since the integrals may not be actually evaluated without 
choosing special velocity distributions, the following 
function is introduced 



Furthermore, corresponding to (35) the velocity integral may 
be calculated as a sum if ia^ v (0) is chosen as a fine-grained 
function. 

Then, by integration of (62) 



Where Z^, z_. , v^, v g ^ are position and velocity components 
of ions and electrons parallel to k at t = 0. 

Equation (64) may be considered the basic equation of 
the theory of longitudinal plasma waves. From this equation 
follow the formulas for scattering, dispersion and drift 
instabilities. 

The next step is to employ the theory of density 

fluctuations to light scattering in a Plasma. This requires 

2 2 

the amplitude squares | Q^(OJ) | and | Q S ^(W) | ; these 

represent a time average and a time integral over all micro- 
scopic states assumed by the plasma between t = 0 and t = 




(63) 



exp ( -ikZ^ ) exp (-ikz.. ) 




-U)-kv. .+iy -u>-k . . +iy 



(64) 
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Substituting (64) into (59) and (61) and integrating 



over all velocities using the same method by which Q] c (^) was 
obtained, | Q s ^(60)| may calculated. 



^ek^ 1-G (05)-G.(o) 



and 

Q ek (t0)= 1-G ^03)- G. (0?) t 



p (1-G i (0)))^exp(-ikz .) £ exp (-ikz.) 

L -a>-kv+iy -40-kv+iy J 



£0-kv+iy 
(65) 



( l“G a ( 03) ) £ exp(-ikZ .) £ exp (-ikz .) 

-co-kv+iy ®i^ ^ -a>-kv+iyJ 

( 66 ) 

The amplitude squares | Q s ^. ( 60) | contain double sums of the 
form. 



£ exp ( ikz . 
j 1 



-o>-kv , +iy 
eD 



£ £ exp [ik (z .-z . ) ] 
= J & £ 1 



(-co-kv e j+iy) (-0)-kv ej2 +iy) 



(67) 



as well as double products. Ion velocities and positions 
are indexed by l, whereas electron velocities and positions 
are indexed by j. Diagonal terms of the left (the single 
sum) are treated the same as off-diagonal terms. In the 
double sum, however, the diagonal terms behave differently; 
that is, the space dependent phase factor vanishes in these 
terms allowing the replacement of the sum with an integral 
over the velocity distribution; 

N 



exp[ik(Z_.-z_. ) ] 
j=l (-0>-kv)^+ y^ ~ ^ 



" F (v) dv 
e ' 

_ co (OH-kv) ^+y 2 



( 68 ) 



Since the denominator of the right hand integral depends 

only on v , the integral is over the one-dimensional 
z 
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distribution F (v) , gotten by integrating F (v) over v and 
Vy. Rigorous evaluation of the integral is at best tedious. 
It may, however, be approximated more readily. The main 
contribution to the integral comes at resonance v = - ^ . 

At resonance, the value of the denominator varies rapidly 
while the numerator varies slowly, permitting it to come 
out of the integral. 



r F (v) dv 

lim j T o 

■y-+0 (a.'+kv) y 



F ( -rr } 
e' k 



dv 



(wt-]cv) 2 + y 2 



’Zefl2 

y k 



( 69 ) 



yielding an expression for the diagonal sum dependent on 
the time of observation 77 . This, of course, is what is 

y 



expected, since | Q^Ca. 1 ) 
tion during T = ^ . 



determines the scattered radia- 



The off-diagonal terms behave somewhat differently. It will 

be shown that this sum may be neglected compared with the 

diagonal terms. The sum is independent of y and may be 

neglected in comparison to the double sum (which depends 

on ~^) for small y. Consider the double sum over a factor 

exp[ik (z .-z . ) ] and assume the denominator is replaced by an 
1 3 

2 

appropriate mean value. The result is a sum over N(N- 1 )-N 
positive and negative terms. From the same "random phase" 
the sum is of the order of magnitude N. Hence, 



II 

j W 3 



exp[ik (z^-z , ) ] 



(-co-kv .+iy) (-u)-kv fl -iy) 
^ e j e l 



- N 



F e ( v n) F e ( v jfc) dv_;dv 



j l 



, (-u>-kv . +iy) (-UJ-kv.-iy) 
3 z 

( 70 ) 
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The single integrals over dv^ and dv^ have linear resonances, 
which by symmetry are integrable; the numerator remains un- 
changed as the denominator changes sign. Thus the resonance 
contribution to the integral is zero. Therefore the single 

integrals are independent of y and may be neglected in com- 

1 



parison to the diagonal terms proportional to 
obtains , 



One then 



Q ek (^ 



=1 n^. ( u>) I = 



I 1-G. (60)1 F hr)+|G (t)rp,w 
ffN i v/l e v k ' e v ' i v k ' 

yk | 1-G e (a')-G i (05)| 2 



( 71 ) 



Q ik (^! 2 



TO I I - G e (*)l 2 F i( ^) + |G i (‘.)| 2 F e r-f) 
yk ! l-G e (0J)-G i (0))| 2 



( 72 ) 



Q je (^) 



F (— ) + F (— ) 

ffN e V k [ l 1 k ; 

7k I l~G e ( 60 ) -G i (U) I 



( 73 ) 



Evaluation of these amplitude squares of course depends on 
calculating G g and G^ for special electron and ion velocity 
distributions . 

Of interest is the shifted Maxwellian distribution; 



F s (v) 



C^kf~) 4 ex P [" W* <2 + w s ) 2 ] 



( 74 ) 
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where W is the drift velocity of particles of type S» 
s 

(S = i,e) 

Clearly, 



dF (v) 

s — 

dv 



r m 1 

[h < v + V J F s (i) 



( 75 ) 



and 



G s<"> 





(v+W g ) F g ( v) dv 
(- co-kv+iy) 



(76) 



where a is given by (30) and G (co) by (63) , and where 
s s 

F (v) is a one-dimensional velocity distribution in v , and 
s z 

W is the z component of W . 
s s 

Introducing the Doppler shifted frequency, one has; 

CO = CO- kW (77) 

sw s 



Quantities of the unshifted (no drift velocity) distribu- 
tion will now be superscripted with a zero. 

Then, 



G s (<4 = 




(v+W )F°(v+W ) d { V+W ) 

s s s s 0 / \ 

» [— co -k(v+W ) +iy] G s W sw 

sw ' s 



(78) 



and 



F (v) dv 
s v 

2 2 

- 00 ( cofkv) +y 




F°(v+W s )d(v+W s ) 
[co +k (v+w ] 2 +y 2 

SW ' s 



JL F o ( 
yk s v k J 

(79) 
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Using the Doppler shifted frequencies 



. ,2 2,22 

to = — k v 
s , s 3 s, s 



2to" 



. s, s ' = i, e 



( 80 ) 



a 



where v , is the mean thermal velocity square of the par- 
s , s 

tides of mass m and temperature T , and introducing the 

S S 

dimensionless frequencies 



x 1 



05 

ew 

05 

ee 



05. 

1W 

05. . 
11 



(81) 



it follows that (using (74)), 



2 2 
l-G i (y 1 ) I exp ( -x 1 2 ) + (G e (x ' ) | T7exp(-y'2) 



ee 



1-G e (x') - G.(y')| 



(82) 



and. 



y* 



05 



05. . 
11 




w. 

1 



v . . 



ii 



U w 

, = 05 _ ( 3' 2 \ 

X 05. . n 05 



11 



ee 



(83) 



where 



V = 




Furthermore, if the substitution 
u = W + v 



(84) 



(85) 
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is made, and the integrand of (79) is expanded with the 
factor (-co-ku-iy) , then. 



G (x) = ka 



2 ’ uF^(u) (-co-ku-iy) du 

-« (-o>ku) 2 +y 2 



( 86 ) 



The imaginary part is an integral of the form of (69) . 
Therefore, 



Assuming y small the real part of G(x) may be expanded in a 
power series, which, when evaluated, is equal to zero in the 
limit. Using the same arguments as with (70) , the real part 
goes through zero at resonance and changes sign. Again, 
integrated over symmetric limits, the integral goes to zero. 

From F^ (u) , let <| u|> be- the mean value of I u|; expand- 
ing the integrand of (86) for CO > k <| u|> in the form 



x 2 

lim I (x) - ax exp(-x ) 

y-0 



(87) 
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-CO-ku CO 



1 . -1 



£(- 1) V <*%) V ( 88 ) 



V=0 



the following is obtained? 



lim Re(x) 
y-*Q 




(89) 



For C0<k <(| u | y , the substitution 



4 = ku - co , u 



U+o>) 

k 



du 



d4 



(90) 



k 
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is made and F (u) is expanded. That is. 



rW(^)-0(-l> xp( - x 2, £=£ (-£-)«* 

V = 0 * KT 



( 91 ) 



Using this to obtain 



lim Re (x) = - a 2 ri-2x 2 (l - | x 2 + Tc x 2 + )1 = 

y-*0 L o Id J 

- Q! 2 Cl-F (x) ] 



( 92 ) 



where f (x) is identified with 



f (x) = Ox exp(-x 2 ) J exp(t 2 )dt 



( 93 ) 



because this function has the same series expansion. Thus, 
finally, 

G s ( x s ') = - a s 2 [l-f (Xg) - iff* x g ' exp ( -x 2 , ) ] s=i, e x/ = y' 

( 94 ) 

What has been presented here is the linear theory for 
scattering. Appendix B gives calculated example curves of 
spectral distributions of plasmas with drift velocity W 
(W is smaller than W^ r , the critical velocity at which non- 
linear drift instabilities are excited) W is in multiples 

3 ^ 

of the electron thermal velocity (tt) 2 v ; a is a , the 

z ee e 
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Each 



T 

parameter identified in equation (30) and e = 7 ^ . 

5 7T 2 ff 

spectrum contains 11 curves from angles <t> = 0 , Yq ' Yo"' 

between W and k with = | W |cos d . The curves with the most 

IT 

pronounced side maxima correspond to <t> = , the only symme- 

TT 

trie curve corresponding to 0 = with the horizontal tan- 
gent at y = 0 . y is of course the dimensionless frequency 

, and extends over | y| =£ 3.5 representing the ion line, 
ii 

Relative maximum intensity I m is given for each, since inten- 
sities vary strongly. Intensities are normalized so that 



I(W,k) 



! 1-G i (y) I 2 exp ( -x 1 2 ) + | G e (x ’)| 2 X] exp(-y 2 ) 
| 1-G e (x') - G.(y) | 2 



(95) 



5 

From the report by Theimer; by permission of 
Professor F. Schwirzke. 
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III. APPARATUS 



The experimental setup itself requires a myriad of 
equipment. The schematic of figure 1 shows the equipment 
involved and its layout. The steady-state plasma facility 
at the Naval Postgraduate School is ideal for a systematic 
investigation of the laser light scattering plasma diagnos- 
tic technique. The parameters of plasma density and elec- 
tron temperature may be varied over wide range in the sys- 
tem. The plasma, whether argon, nitrogen, or helium, is 
produced through a hollow tantalum cathode and is acceler- 
ated to a disc anode between voltages of around 40 volts. 

The beam is columnated at the cathode by an arrangement of 
wound electromagnets producing fields of several hundred 
gauss. The plasma density may be controlled by varying the 
gas flow at a needle valve, altering the background pres- 
sure with a 10 inch diffusion pump, and by adjusting the 
field strengths of the columnating electromagnets. Plasma, 
and hence electron, densities^ asj, seen in the theory, have 
a direct bearing on the scattered light intensities. Pre- 
paration of the system required the construction and place- 
ment of perpendicular window ports through which the inci- 
dent light entered and the scattered light emerged. The 
ports were approximately a half inch in diameter and were 
mounted on removable flanges, as shown in figure 6. The 
plane described by the perpendicular incident and scattered 
light wave vectors was arranged in such a way as to intersect 
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the plasma jet at a point within a half inch of the tanta- 
lum cathode outlet, for two reasons. One, this is presum- 
ably where the beam and hence the electron density are the 
greatest; secondly, little is known about the plasma para- 
meters at the cathode... the scattering technique should 
yield much information concerning these parameters. Further- 
more, laboratory space limitations and apparatus configura- 
tion demanded location at the spot chosen. 

The incident laser radiation was provided by the model 
130 high peak power Q-switched laser, manufactured by Optics 
Technology, Inc. The instrument consisted of a power 
supply/control console and a separate optical head, which 
incorporated a Fabry-Perot type interferometer and laser 
cavity mounted on a granite base. The Q-switch is a device 
which enables the laser to deliver extremely high peak powe'r 
pulses varying from one to a hundred megawatts. As here em- 
ployed, it was a bleachable liquid system contained in a 
small cell with optically flat faces, which only needed to 
be inserted between the laser ruby and the interferometer 
for operation. Further function and use of this device 
will be discussed later. The laser crystal is a production' 

quality, 60° orientation, three-eighth by three inch. ruby 

o 

crystal, emitting radiation at 6943 A. The crystal is op- 
tically pumped by two linear xenon arc lamps with three 
inch arc spacing and overall diameter of eleven mm, mounted 
opposite the crystal holders in a polished aluminum cavity. 
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Cooling is accomplished through a forced draft, water-cool- 
ed air system. The rear element of the Fabry-Perot inter- 
ferometer is a TIR prism, and a multi-element glass re- 
flector mounted on a micrometer driven assembly comprises 
the front element. Alignment of the interferometer optics 
is extremely critical for lasing action to occur. Align- 
ment may be accomplished in the following way. With the 
top of the laser cavity and the ruby crystal removed, care- 
fully position a cw laser so that its beam is in line with 
the axis of the ruby crystal holder pins and passing through 
the center of the exit hole as shown. 



RUBY LASER 



ruby holders 
flash lamps 



adjustable / 
mount 



CW LASER 

•--] = 



The task of aligning the cw laser may be more easily accom- 
plished accurately by placing a pen or a pencil on the ruby 
pins, holding it vertically, and adjusting the cw laser 
until the beam strikes the pencil in the center when placed 
on either pin. There will now be a multitude of reflec- 
tions on the face of the cw laser, one from the TIR prism 
and one from each of the glass reflectors in the microme- 
ter driven interferometer. Determine which one of these 
reflections comes from the ITR prism by placing a sheet of 
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paper directly in front of the prism; Once finding the 
correct reflection, remove the paper and hand adjust the 
prism until its reflection is as close to the cw exit aper- 
ture as possible. Then tighten its screws at its base to 
secure it. Next adjust the micrometer screws on the inter- 
ferometer until the multiple reflections made by its ele- 
ments merge on or near the cw exit aperture and prism re- 
flection. This accomplished, fringe patterns should be 
observed, and the ruby laser is aligned. 

The bleachable dye Q-switch device previously mention- 
ed permits the formation of giant pulses. The technique, 
usually referred to as passive Q-switching employs a small 
cell containing a solution of one of several metallo-organ- 
ic compounds known as phthalocyanines . In this case, a 
solution of cryptocyanine in methanol is used. The- solu- 
tion strongly absorbs light of the ruby frequency, and this 

absorption prevents net amplification of light until a 

3 + 

greater proportion of Cr ions have been pumped into the 
upper energy state than is normal. The energy input con- 
tinues until the amplification in the ruby overcomes losses 
due to absorption in the cell. To this point the laser 
smits coherent light weakly. A small amount of the light 
bleaches the cryptocyanine which then becomes nearly per- 
fectly transparent to the ruby light. At the same instant 
there is a large net amplification and a giant pulse con- 
taining all the stored ruby energy develops rapidly. 
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Following the pulse, the solution returns quickly to its 
absorbing state. Below is a typical phthalocyanine absorp 
tion curve: 




Bleaching of the metal cryptocyanines is effected by the 
saturation of the absroption at the ruby frequency. That 
is, the solution is in the bleached condition whrn the pop 
ulation densities in the upper and lower states of the 
phthalocyanine molecule (S',S in the diagram below) are 
roughly equal. Because of the large cross-sections for 
capturing ruby light only small numbers of cryptocyanine 
molecules are involved. Consequently saturation occurs 
rapidly and relatively few quanta from the laser beam are 
required. 
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excited t riplet T ’ 




The solution itself is prepared from the concentrate 
by diluting one part cryptocyanine to ten parts methanol. 
The cell is filled almost to capacity with methanol and 
nine drops of the 1:10 solution is added to this. This 
appears to give optimum Q-switching without requiring too 
high a threshold. 

Calibration of the laser output, that is, determining 

output energies and pulse widths (and hence output power) 

involves the use of a PIN photodiode (See Appendix C) and 

a laser radiometer. The HPA 4205 is a silicon planar PIN 

photodiode with a response time of less than one nanosecond 

9 

and a frequency response from dc to 10 hz. The diode was 
wired according to the following diagram. 
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PHOTODIODE 



20 V ^=1 



The photodiode is useful in conjunction with a fast-rise 
oscilloscope, such as a tektronix 519, for recording laser 
pulse shapes. 

Laser output energy was measured using a Westinghouse 
model RN-1 Radiometer which is, more specifically a calor- 
imetric bolometer. It is designed to trap the energy beam 
from the laser in a bundle of fine, insulated copper wire 
which functions both as a calorimetric m ass and as a bolo- 
metric element. The change in resistance of the wire ele- 
ment is proportional to the energy absorbed and is almost 
entirely independent of the distribution of the energy with 
in the unit. The bolometer element is part of a convention 
al Wheatstone bridge circuit and the change in balance of 
the bridge is a measure of beam energy. (See figure 4.) 
When used in conjunction with an electronic microvoltmeter, 
laser energy may be read directly. When a microvoltmeter 
reads the imbalance in microvolts, the reading is directly 
proportioned to the incident laser pulse energy. 




in microvolts 




^ = Laser output in joules 

(94) 
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The schematic of the calibration setup is as follows: 




osci I lose ope 



The following are photographs taken on a Tektronix 
545A oscilloscope set to trigger as the pulse was initiated 
The time scale on the three non Q-switched photos was 504 
s/div. , and on the fourth, Q switched pulse 500 ns/div. or 
0.54s/div. This gives an idea of the effectiveness of in- 
creasing lasing power by greatly shortening the pulse 
width. The Q-switched pulse width appears to be of the or- 
der of 200ns. The average pulse length of the non Q-switch 
ed pulse appears to be of the order of 2504s. Since power 
is energy per unit time, this shows an increase in lasing 
power by a factor of 1250. More specifically, a typical 
low power pulse will put out approximately 0.5 joules. For 
a pulse width of 2504s this comes out to two kilowatts. 

The Q-switched pulse width of 0.24sec. the power is now 
2.5 megawatts, a substantial increase. The first three 
photographs show 3 pulses each; below them is recorded 
the average energy output of the three for a given input 
energy from the laser power pack. This was the same for 
each trace on a given photograph. 
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Non Q-Switched Pulses 
Sweep Speed = 50 Jls/div. 




Input 1100 joules 
Output 0.19 joules 



Input 1225 joules 
Output 0.38 joules 
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Non-Q-Swi tched 
Sweep Speed 50 jis/div. 




Input 1300 joules 
Output 0.55 joules 



For these first 3 photos, a plot of input versus output 
energy is made, establishing a linear relationship above 
laser threshold. 




1100 1200 1300 1400 

Power Pack Energy Input (Joules) 
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The fourth photo is included as a fine example of a 
Q-switched laser pulse. Sweeptime here is 0.5 (ls/di vision. 




The scattered light is analyzed with the Baird atomic 
modified Eagle Mount three meter concave Grating Spectro- 
graph. The spectrograph is a Rowland circle type mounting 
for which the focusing conditions are: 



r = RcosA A 

r'= RcosB B 



( i 




= Angle of incidence 
= Angle of diffraction 
= Slit to grating distance 
= Camera to grating 
distance 

= Radius of curvature of 
the concave grating. 

= Normal 
= Grating 
= Slit 
= Camera 
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The Rowland circle has a diameter equal to the radius 
of curvature of the grating blank. The grating is mounted 
tangent to the Rowland circle. 

The usual grating equation applies. That is 

nX - d(sin A + sin B) (95) 

n = Order number 
X = Wavelength in Angstroms 
d = Grating spacing in Angstroms 
A and B are of opposite sign if they are on different 
sides of the normal. Spectrograph controls include toggle 
switches for the grating, the plateholder (on which the 
photomultiplier tubes are mounted) and focus. 

Rotating the grating changes the angle of incidence 
A (and hence B) . The Veeder-Rcot counter is geared so as 
to indicate the angle of incidence, subtract from 1000. 

The focus toggle increases or decreases the plateholder to 
grating distance, the tilt control moves the plate (mount- 
ing the photomultipliers) in and out about a vertical axis 
through the center of the slit. The latter two controls 
are not as important when using the photomultiplier as op- 
posed to a photographic plate when a careful focus is 
required. 

There existed no calibration chart of the instrument 
to predetermine the position of the photomultipliers for a 
specific wavelength. This must be accomplished experiment- 
ally, using a calibration source of known wavelength, for 
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example a helium lamp* Helium has two lines close to the 
6943 A ruby wavelength, one at 6730 A and another at 7108 A. 
Sighting visually through the plate aperture, these may be 
readily identified. With a slit cut in dark cardboard, the 
position at which these lines strike the aperture may be 
marked. Linearly interpolating between the two wavelengths 
gives a rough position for placement of the photo multi- 
plier. A fine adjustment may be made when the whole system 
is aligned. The alignment procedure will be described 
shortly. For purposes of this experiment the proper grat- . 
ing setting was found to be 991.2 (that is, with the angle 
of incidence A equal to 8.8 degrees.). The dispersion may 
also be calculated with some ease. This will become of 
great use when actual spectral distributions are analyzed 
carefully. This relation may be derived from the basic 
grating equation. 

nX = d(sin A 4 sin B) (95) 



X = —(sin A + sin P) 

n 



letting d be a constant 




B 



(96) 



or with A constant 



dX 

dB 



— cos B 
n 
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(97) 



an = 



dB 
d l 



dX 

d£ 



dX 

d£ 



but 

r = 



hence 



d X 
d£ 



dX 

d l 

d = 

therefore 

dX 

dJL 



rdB 

where r is the focal distance and 
is the increment along circumference 
for dB increment. 



(98) 

dX . dB 
dB * an 



- cos B 
n 



1 

r 



(99) 



R cos B (100) 



d cos B D = actual grating spacing. 

nR cos B 

d = effective grating spacing. 



d. 

nR 



( 101 ) 



D cos B 



( 102 ) 



D cos B 
nR 



(103) 
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Numerical calculations for dispersion may now be made, let 

ting D be the spacing for 15000 groove/inch grating and 

o 

R = 3 meters at the 694 3A~ wavelength and an angle of inci 
dence A (determined experimentally) of 8.8°. 



D 



1 

15000 



inches 



254mm 10 A 

inch x cm 



16933A 



for the first order spectra 



X = d(sin A + sin B) 



A = 8.8 



sin B = — - sin A 
d 



6943 . „ O o 

16933 ' Sin 8 ’ 8 



= .41003 - .15299 

= .25704 
cos B = .96638 



dX _ D cos B 

an r- 



(16933) (.96638) 
3 x 10 3 



16369 
3 x X0 3 



o 

5.45 A/mm 
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The photomultiplier tubes were RCA 7102 tubes with 



sensitivity and response curves given below. 




wavelength 
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sensitivity amperes/ lumen 




supply volts 



The tubes were operated from a number of power sup- 
plies, the best of which appeared to be the Power Designs 
Model HV 1565 regulated dc power supply. Optimum signal 
to boise ratio was obtained at a voltage setting between 
1450-1490 volts (actual voltage as read on the voltmeter 
was closer to 1350 volts) . The output signal was read 
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current amplification 



across about a five megohm impedance for greater signal 
strength. This was of course at some sacrifice to response 
time, as will be seen later. 

Before discussing the important alignment procedure, 
the optics used to transfer the scattered light from the 
plasma to the spectrograph will be described. The arrange- 
ment was a periscope type as shown in figures 1 and 5. Two 
converging lenses were employed, one focusing at the plasma, 
the other refocusing the plasma light at the exit of the 
periscope. A quartz condensing lens on the optical bench 
of the spectrograph refocused the light on the slit. An 
arrangement of rigid cardboard packing tubes and aluminum 
mirror holders formed the remainder of the device. The po- 
sition of the periscope exit (on the plane of which the 
plasma beam is focused) is determined by using the grating 
lamps and quartz condensing lens in the following way. The 
grating lamps are hung over the front face of the grating 
holder and switched on from a control on the external power 
supply. The quartz condensing lens is placed in its holder 
on the optical bench and left free to move along the length 
of the bench. The image of the grating lamps may be observ- 
ed by holding a white sheet of paper behind the lens. Mov- 
ing the paper away from the lens will eventually focus the 
grating lamps in two distinct images. Moving the paper 
further away, the images will merge, and it is at that 
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point that the source, in this case the exit of the peri- 
scope, is to be placed. The lense is now secured and the 
periscope is mounted in its proper position. 

Because of the distance it was necessary to convey 
the lights and because of the small image of the lenses, 
alignment was extremely critical and somewhat difficult to 
obtain. Accurate laser-plasma-periscope-spectrograph align 
ment may be obtained while the plasma system is at atmos- 
pheric pressure and the tantalum cathodes are being replaced 
With the cathode flange removed, remove also the left win- 
dow flange. Place a wooden slate across two stationary 
bolts of the cathode flange in such a way as to have the 
slat horizontal and just below the center of the vacuum 
chamber (and hence below the plasma beam axis) . The care- 
fully insert a rigid aluminum or stainless steel rod down 
the chamber, supporting one end inside the annular disk 
anode and the other end on the wooden slat. The slat 
should be of such thickness that the axis of the rod it 
supports is parallel to and just below the chamber axis. 
This arrangement will now support the mirror proble to be 
inserted from the left window opening; the probe will act 
as an artificial "scatterer" so that proper alignment both 
of the optical transfer system and of the grating may be ob 
tained. The mirror probe consists of a straight wooden 
dowel with one end cut to 45 > . To the 45° cut is applied 
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a small cutout piece of polished stainless steel to act as 
a reflector. This may be glued on and a mark may then be 
made on the opposite end of the dowel to indicate the pro- 
per vertical mirror position. The dowel is then inserted, 
one end supported on the rod arrangement described above 
and the other end on a notched cardboard edge, taped to 
the window opening. The supporting rod may vary the position 
of the reflector along the chamber axis without significantly 
changing its position perpendicular to the axis; pushing 
or pulling the mirror rod itself, varies its position per- 
pendicular to the chamber axis . 

The next step is to place a strong light source such 
as a mercury lamp before . the grating (inside the specto- 
graph) and aimed back toward the slit. The Mercury beam may- 
be focused at the slit using its own lens for more intensi- 
ty. With the laboratory darkened, carefully adjust the per- 
iscope arrangement until the light from the Mercury falls 
on the supporting rod, which has been located parallel to 
the vacuum chamber center axis. This done, the mirror prop 
may be adjusted so that the light falls on its reflector. 

From the end of the mirror proble, and peering into the 
chamber through the left window mount the probe may be ad- 
justed carefully so that its light falls on the half inch 
laser window. The axis of the mirror should be vertical, 
that is, reflecting straight up. With the Mercury light, 
(which may be easily discerned) shining on the half inch 
laser entrance window, the system, save the laser, should 
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now be in perfect alignment. The next problem is to align 
the laser so that its artificially "scattered" light may be 
seen by the grating; and when this is done the grating must 
be rotated until the laser frequency is centered on photo- 
multiplier tubes at the spectrograph exit. 

The laser is mounted on an aluminum vertically adjust- 
able stand, the laser itself containing tilt adjustment on 
its base. Alignment along the horizontal plane may be made 
more easily with a cw laser by raising the axis of the cw 
laser, the pulse laser may replace along the same axis and 
trial and error adjustment made until alignment in the ver- 
tical plane is achieved at will. Once aligned the laser 
pulse should be seen in the slit when fired. 

The grating may then be rotated to the proper wave- 
length. The phototubes and grating have already been rough- 
ly adjusted according to the helium line procedure described 
above. First the PIN photodiode is positioned with an arran- 
gement of beamsplitters (microscope slides) to analyze the 
laser pulse. The signal from the PIN photodiode is then fed 
to the upper trace of a dual beam oscilloscope, and the sig- 
nal from the photomultiplier to the lower trace. Upper 
trace gain is lowered because of the intensity of the laser 
pulse, but the oscilloscope arrangement is as shown in 
figure 14. The grating is rotated by 0.1 degree at a time 
in order to find the optimum position. Below are the traces 
illustrating how this was accomplished. The upper trace is 
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the pm tube signal and the lower trace is the laser pulse. 
Beside each is the respective grating setting and its cor- 
responding angle of incidence. 




991.2 (8.8°) 
991.1 (8.9°) 
991.0 (9.0°) 
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991.5 (8.8°) 



991.3 (8.7°) 



991.2 (8.8°) 



Sweepspeed was 50 jiis/div. Both traces were triggered by 
the laser pulse, the lower trace of each group. The laser 
trace gain was 0.5 V/div. and the pm tube trace 25 V/div. 
All were non Q-switched pulses. Note the difference in 
pulse shapes. As was mentioned before, this was due to 
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the large resistor in the signal line (about five megohms) 
which considerably lengthened the rise time of the tubes. 
This of course, is no consequence since relative intensi- 
ties are what is important here and in analyzing the spec- 
tral distribution of the scattered light. 

The system, now constructed and aligned, may be em- 
ployed in the attempt to discern a scattered light signal 
and analyze the spectral distribution of the scattered 
light. The few attempts thus far have been unsuccessful. 
Some combinations of detecting tubes and oscilloscope am- 
plifiers have been tried to increase the signal to noise 
ratio. Photographs have been taken, an example of which 
is below. 
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The photograph on the bottom was triggered at the 
scope with no laser pulse in order to provide a noise refer 
ence from the pm tube. Tube voltage was set at 1490 (Ac- 
tually 1390 on the voltmeter) . Sweep speed was 50 /ls/div. 
and the background pressure in the vacuum chamber was 
8 x 10 ^ mm Hg. Oscilloscope gain here was 0.05 volts/div. 
as shown in figure 14. If a signal is present, as a few of 
the traces such as the one above seem to indicate, it is 
virutally impossible to distinguish it from the noise. 
Possible improvements will be discussed in the conclusion. 

Attempts at focusing the laser beam seem to have no 
effect, but this is almost certainly due to the poor quali- 
ty of lenses available. 
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CONCLUSIONS 



IV. 

The scattering apparatus has been erected and aligned. 
The attempts which have thus far been made, to discern the 
scattered light and analyze its spectrum have been unsuccess- 
ful. The plasma parameters specifically the density, affect 
scattering intensities as much as many varied instrument 
parameters affect the experimenter's ability to detect the 
light. All oscilloscopes, save a very few are in dire need 
of realignment and recalibration, as are the plug-in units 
that accompany them. When the scattered light is finally 
detected up to seven photomultipliers will be needed to si- 
multaneously analyze the spectral distribution of a single 
pulse; knowing the dispersion of the grating, the photo- 
multipliers, connected to the spectrograph plate with short 
lengths pipe, may be accurately spaced around the central 
frequency to observe the whole spectral distribution of the 
scattered pulse. These should be of the RCA 7102 type. 

The 7102 tubes now available appear to have lower signal 
to noise ratios at sensitivities required that one would 
expect. Improvement of that ratio could be of significant 
importance in detecting the scattered light. Increasing the 
plasma density by purging more gas into the plasma system, 
increasing the field on the cathode and anode magnets and 
by controlling the main diffusion pump may well be the 
answer sought. Furthermore, the laser system available. 
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while powerful, may not be quite what is needed, though 

focusing of its beam through a lens of better quality may 

help. A more energetic system is now on order. It is 

interesting to note previous experimental attempts to study 

at light scattering and the problems mentioned in these 

efforts. One group, Thompson and Fiocco at MIT measure the 

light of a laser scattered from an electron beam. The elec- 

9 -3 

tron density was about 5 x 10 cm , and the ratio of 

— 18 

scattered to incident intensity was of the order of 10 1 

Applying their technique to a low electron density plasma, 

13 -3 

n « 10 cm , they improved the ratio of scattered to in- 

-15 

cident power to 10 

The best scattered to incident ratio came from E. 

Funfer, B. Kronas t, and H. J. Kunze, who measured light 

17 

scattered from a high density pinch plasma (n « 10 ) . 

-11 

Their ratio came to 10 

Distinguishing the scattered light from background 
plasma radiation is a considerable problem. One possible 
solution that may be tried is to eliminate the plasma ra- 
diation with polarizers, since at 90° the scattered light 
will be polarized, whereas the plasma background radiation 
will not. One other wasy of eliminating plasma radiation 
would be with the use of a differential amplifier which will 
differentiate the scattered signal from a pair of pm tubes 
one seeing the scattered light and radiation, the other only 
radiation. This technique was attempted, but the character- 
istics of the two pm tubes available were not close enough 
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for comparison. Furthermore, the differential amplifier 
plug-in units on the oscilloscopes were not functioning 
properly. 

The above has summarized experimental progress to 
date. This progress should enable future experimenters to 
resume development of the experimental technique of laser 
light scattering as a plasma diagnostic. 
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APPENDIX A. DERIVATION OF THE BOLTZMANN EQUATION 6 

Each electron is represented by a point in six dimen- 
sional electron phase space. The coordinates of the phase 
space are x, y, z, u^, u , u z . Because each electron has 
a velocity, the representative point for each electron must 
change its position coordinates with time. In addition, 
any change in the velocity of the electron will also cause 
the representative point to move in electron phase space. 
The velocity of an electron will change if a force acts on 
the electron or if the electron collides with any other 
particle. 

Consider a fixed differential volume in electron phase 
space dV = dx dy dz du du^ du z (the "electron" subscript 
will be omitted) . An expression will be written for the 
number of representative points in this volume as a func- 
tion of r_, u, and _t; that is, it will be a "continuity" 
equation for the representative points in electron phase 

space. The increase with time in the number of electron 

3(f e dV) 

"points" in dV, §■- is caused by three factors: 

1. The drift of points in position coordinates due 
to inherent velocity. 

2. The drift of points in velocity coordinates due 
to electric and magnetic forces. 



6 

From Uman, p„ 34-36. 
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3. The spontaneous velocity change without position 
coordinate change due to collisions. 

When these three terms are evaluated, the Boltzmann equa- 
tion for electrons may be written. Note that the effects 
of external gravitational forces and inteparticle nuclear 
forces are omitted. 

First determined is the increase with time of the 
number of representative points in dV due to drift in po- 
sition coordinate space. The number of representative 
points which cross the surfaces of dV in time de because 
of inherent velocity must be calculated. The velocity of 

a point in the x direction at x, y, z, u, u , u is u . 

y z x 

The number of representative points in a volume of "base 

area" dy dz du du du and height u dt will enter dV in 

time de through a face normal to the x axis at x, y, z, 

u , u , u . This number is 
x' y z 

[ ] dt = f Q (r,u, t)u x dtdydzdu x du du z (1) 

X 



The change in the number of representative points in dV in 
time dt due to the exit of points through the equivalent 
surface at x+dx is given by 



3(d 6 N) ~ 

-* x+dx 



dt = 



-[ 



3 (f u ) 



f u + 
e x 



e x 



3x 



dxj dt dydzdu^du du r 



( 2 ) 
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